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Sho Matsumoto and Tomoyuki Shirai 

Abstract 

We show that the zeros of the random power series with i.i.d. real Gaussian 
coefficients form a Pfaffian point process. We also show that the product moments 
for absolute values and signatures of the power series can also be expressed by 
Pfaffians. 



1 Introduction 

Zeros of Gaussian processes have attracted much attention for many years both from 
theoretical and practical points of view. The first significant contribution to this study 
was made by Paley- Wiener [18]. They computed the expectation of the number of zeros 
of (translation invariant) analytic Gaussian processes on a strip in the complex plane that 
are defined as Wiener integrals. Their work was motivated by papers due to Bohr-Jessen 
[1, 2] on almost periodic functions in the complex domain arising from Riemann's zeta 
function. Kac gives an explicit expression of the probability density function of real zeros 
of a random polynomial 



k=0 

with i.i.d. real standard Gaussian coefficients {a^j^^Q and obtains a precise asymptotics 
of the numbers of real zeros as n — oo [11]. Rice also obtained similar formulas for the 
zeros of random Fourier series with Gaussian coefficients in the theory of filtering [19]. 
Their results have been extended in various ways (e.g. [5, 14, 20]) and generalizations of 
their formulas are sometimes called the Kac-Rice formulas. A recent remarkable result 
on zeros of Gaussian processes is that the complex Gaussian process fc{z) := J2T=o(kZ^ 
with i.i.d. complex standard Gaussian coefficients form a determinantal point process on 
the open unit disk © associated with the Bergman kernel K{z,w) = , which was 

found by Peres- Virag [17]. Krishnapur extended this result to the zeros of determinant of 
the power series with coefficients being i.i.d. Ginibre matrices [13]. 
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In the present paper, we deal with the Gaussian power series 

oo 

f(z) = J2(^kz', (1.1) 

fe=0 

where {ak}'^Q are i.i.d. real standard Gaussian random variables. The radius of conver- 
gence of / is almost surely 1, and the set of the zeros of / forms a point process on the 
open unit disc © as that of fc does. The primary difference between / and fc comes from 
the fact that f{z) is a real Gaussian process when the parameter z is restricted on (—1,1) 
and each realization of f{z) has symmetry with respect to the complex conjugation so 
that there appear both real zeros and complex ones in conjugate pairs. The point process 
of real zeros and that of complex zeros are singular each other. 

Our main purpose is to show that both correlation functions for real zeros and com- 
plex zeros of / are given by Pfaffians, i.e., they form Pfaffian point processes on (—1, 1) 
and D, respectively. The most known examples of Pfaffian point processes appeared as 
random eigenvalues of the Gaussian orthogonal/ symplectic ensembles. Real and complex 
eigenvalues of the real Ginibre ensemble are also proved to be Pfaffian point processes 
on M and C, respectively [3, 7]. Recently, it is shown that the particle positions of in- 
stantly coalescing (or annihilating) Brownian motions on the real line under the maximal 
entrance law form a Pfaffian point process on M [21], which is closely related to the real 
Ginibre ensemble. Our result on correlation functions of zeros of / is added to the list of 
Pfaffian point processes, which is also obtained independently in [6] via random matrix 
theory. Here we will give a direct proof by using Hammerslcy's formula for correlation 
functions of zeros of Gaussian analytic functions and a Pfafhan-Hafnian identity due to 
Ishikawa-Kawamuko-Okada [10]. This is a similar way that was taken in [17] to prove that 
the zeros of fc form a determinantal point process, and in the process of our calculus for 
real zero correlations, we obtain new Pfaffian formulas for a real Gaussian process. The 
family {/(i)}_i<t<i can be regarded as a centered real Gaussian process with covariance 
kernel (1 — st)~^. We show that, for any —1 < ti,t2, ■ ■ ■ ,tn < 1, both the moments of 
absolute values E[\f{ti)f{t2) ■ ■ ■ /(^n)|] and those of signatures E[sgn. f{ti) ■ ■ •sgn(/„)] are 
also given by Pfaffians. We stress that it should be surprising because such combinatorial 
formulas cannot be expected for general centered Gaussian processes. These are special 
features for the Gaussian process with covariance kernel (1 — st)~^. 

The paper is organized as follows. In Section 2, we state our main results for correal- 
tions of real and complex zeros of / (Theorems 1 and 6), and we give new product moment 
formulas for absolute values and signatures of / (Theorems 4 and 5). Also we observe 
negative correlation property of real and complex zeros by showing negative correlation 
inequalities for 2-correlation functions. The asymptotics of the number of real zeros in- 
side intervals growing to (—1, 1) is also shown. In Section 3, we recall the well-known 
Cauchy's determinant formula and the Wick formula for product moments of Gaussian 
random variables. In Section 4, after we show an identity in law for / and /' given that 
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/ is vanishing at some points, we give a preliminary version of Pfaffian formulas (Propo- 
sition 4.4) for the derivation of the expectation of product of sign functions. In Sections 
5, 6 and 7, we give the proofs of our results stated in Section 2. 



2 Results 

2.1 Pfaffians 

Our main results will be described by using Pfaffians. Let us recall the definition of 
Pfaffians. For a 2n x 2n skew symmetric matrix B = {bij)'^'^^i, the Pfaffian of B is defined 
by 

Pf(^) = X](sgn^)&»,(l)r,(2)^(3)7,(4) • ■ ■br,{2n-l)r,{2n), 

summed over all permutations 77 on {1,2, ...,2n} satisfying r]{2i — 1) < r]{2i) {i — 
1,2, . . and 77(1) < 77(3) < • • • < 77(271 — 1). Here sgn77 is the signature of 77. For 
example, 

Pf(S)=&n if 71 = 1 and Pf(-B) = 612634 - &13624 + &14&23 if ti = 2. (2.1) 

For an upper-triangular array A = (aij)i<j<i<2n; we define the Pfaffian of A as that of 
the skew-symmetric matrix B = (6ij)ij=i, each entry of which is bij = —bji = aij if i < j 
and bii — 0. 



2.2 Notation 

For — 1 < s,t < 1, we define the functions 



1-st' ^cr(s,s)a(t,t) 1-st ^ ' 



and the skew symmetric matrix kernel IK by 

K{s,t) = 

with 



Kn(s,t) Ku{s,t) 

K2l(s,t) K22(s,t) 



Kn{s,t)^^= -, Ki2(s,i) ^ ^ 



V(l-s2)(l-t2)(l_5t)2' Vl-Sn-Si' 

jl - 52 I 

IK2i(s,i) = - Y i_^2i_ 1^22(5, i) =sgn(i - s) arcsinc(s, i). 
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(2.3) 



where sgn(t) = \t\/t for t and sgn(t) = for i = 0. Note that Ki2(s, t) = -K2i(i, s) 
and 

^ V|^22(^,t) K22(5,t) )■ 

We denote by Pf (]K(tj, tj))^^-^]^ the Pfaffian of the 2n x 2n skew symmetric matrix 

fK{ti,h) K{h,t2) ... K(ii,i„)\ 
K(i2,ii) K{t2,t2) ... K(i2,y 



\^K(i„,ii) K(t„,t2) 



K(i„,i„)/ 



Thoughtout this paper, j£„ denotes the set of all sequences t = {ti, . . . ,tn) ofn distinct 
real numbers in the interval (—1, 1). 



2.3 Real zero correlations 

Our first theorem states that the real zero distribution of / forms a Pfaffian point process. 

Theorem 1. Let pn be the n-point correlation function of real zeros of the Gaussian power 
series defined in (1.1). Then 

Pnitl, ...,tn)^ TT-" Pf (K(ii, tj))lj^^ (-1 < ti, . . . ,tn < 1). 

From the explicit expression of correlation functions, for example, the first two corre- 
lations are given as follows: 

pi{s) = 7r"^Ki2(s,s) 

P2{s, t) = 7r-'{Ki2(s, s)Ku{t, t) - Kn{s, t)K22{s, t) + K^is, t)K2i{s, t)}, (2.4) 
from which we easily see that 

as t — )■ s. The first correlation is observed by Kac and many others although Kac consid- 
ered the random polynomial with i.i.d. real Gaussian coefficients. The second asymptotic 
expression means that the real zeros of / are repulsive each other as expected. Moreover, 
we can show that the 2-correlation is negatively correlated. 

Corollary 2. Let R(s.t) = be the normalized 2 -point correlation function. Then, 

R{s, s) = 0, R{s, ±1) = 1 and R{s, t) is strictly increasing (resp. decreasing) for t € [s, 1] 
(resp. t G [— l,s]^. In particular, P2{s,t) < pi{s)pi{t) for every s,t E (—1, 1). 
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By using (2.4), we can also compute the mean and variance of the number of points 
inside [— r, r]. 

Corollary 3. Let be the number of real zeros in the interval [— r, r] for < r < 1. 
Then, 

ENr = - log ^i^. Var A^^ = 2 f 1 - ENr + 0(1) 

TT 1 — r ' V TT/ 

05 r — > 1. 

Renicirk 2.1. The kernel K in Theorem 1 is not determined uniquely. For example, we 
can replace K by K', which is defined by 

ttW / N Sgn(t — S) . / N 

K22{s,t) = == axcsm c{s,t). 

In fact, if we set 

fVT^ 

Q[s,t) = dst 







then {Q{t,,tj))lj^^ ■ {K{ti,tj))l^^^ ■ {Q{ti,t,))l^^^ = {K'{ti,tj))l^^^, and therefore two 
Pfaffians associated with IK and K' coincide from the following well-known identity: for 
any 2n x 2n matrix A and 2n x 2n skew symmetric matrix B, Pi{ABA^) = (det ^)(Pf B). 

2.4 PfafRan formulas for a real Gaussian process 

As corollaries of the proof of Theorem 1, we obtain Pfaffian expressions for averages of 
• ■■f{tn)\ and sgn/(ti) • • •sgn/(y . 

Theorem 4. Let f be the Gaussian power series defined in (1.1). For t — {ti, . . . ,tn) £ 
Xn, we have 

E[\m)f{t2)---f{tn)\]^(^l^ ' (detE(t))-^Pf(K(i„t,))^,=„ 

where E(t) is the n x n positive- definite symmetric matrix given by 

1 



m-Hu,tj))i.^,, C7{s,t) = j 



- St 



6 



Matsumoto and Shirai 



Theorem 5. Let f be the Gaussian power series defined in (1.1). For {ti, . . . , t2n) £ ^2n, 
we have 

E[sgn f(ti) sgn /(ta) • • • sgn f(t2n)] 

^) n sgnfe-t.)-Pf(K22(t„t,))2=,. 

' ^ l<i<j<2n 

In particular, if —1 < ti < t^ < ■ ■ ■ < t2n < ^ , then 

E[sgnf{t^)sgnf{t2)■■■sgnf{t2n)]= 0^ Pf (arcsin c{ti, ^j))i<j<j<2„ > 

(2.5) 

where c{s,t) is defined in (2.2). 

We note that it is easy to observe E[sgRf{ti) ■ ■ ■sgn/(t„)] = when n is odd. More 
generally, if (Xi, . . . , X„) is a centered real Gaussian vector, then E[sgn Xi ■ ■ ■ sgn X„] = 
for n odd. Moreover, the formula (2.5) can be rewritten as 

E[sgnfit,) sgn /(t2) ■ ■ ■ sgn /(t2„)] = Pf {E[sgnf{t,) sgn /(t,)])i<,<,.<2„ • 

As explained later, the Wick formula (3.6) provides us a similar formula for products of 
real Gaussian random variables, however, such neat formulas for £'[|XiX2 • • • are not 
known for general n except the cases with n = 2,3([15, 16]). Similarly, there is no known 
formula for E'fsgn Xi sgn • • • sgn X2n] except the n — 2 case 

£'[sgn Xi sgn X2] — — arcsin ■ '^^^ 



TT \J(y\\(y22 ' 

where = E[XiXj] for i,j — 1,2. Theorem 4 and Theorem 5 state that the moments 
E[\Xi - ■ ■ Xn\] and i?[sgnXi ■ ■ ■ sgnX„] have Pfaffian expressions if the covariance matrix 
of the real Gaussian vector (Xi, . . . , X^) is of the form ((1 — titj)~^)2j^i. 

2.5 Complex zero correlations 

The complex zero distribution also forms a Pfaffian point process. Put D+ — {z & 

C I 1^1 < 1, '^z > 0}. 

Theorem 6. Let be the n-point correlation function for complex zeros of f . For 

Zi, Z2-, ■ ■ ■ , Zn & D_|_, 



n 

■ ■ ■ , ^n) = , n rr^T^ • ^^■))^^=i' 

where 'E^{z,w) is 2x2 matrix kernel 

K\z,w) -- 



(1 — 2U>)2 (1 — 2U>)2 
Z—W Z—W 

{l—zw)"^ (1— ztu'l^ 
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For example, the one-point correlation is 



and 



\z 


— z\ 










7r|l-z2 


(1- 


z 








1 




( 


z — W 


2 


z — w 








1 — zw 




1 — zw 



It is easy to verify that P2{z,w) < p'{{z)pl{w) for z,w & D+, which implies negative 
correlation as well as the case of real zeros. 

As we mentioned, Theorem 1 and Theorem 6 are obtained independently in [6] via 
random matrix theory, but Theorem 4 and Theorem 5 are new. 



3 Cauchy's determinants and Wick formula 

In this short section, we review Cauchy's determinants and Wick formula, which are 
essential throughout this paper. 

3.1 Cauchy's determinant and its variations 

The following identity for a determinant, the so-called Cauchy determinant identity, is 
well known in combinatorics. See, e.g., [4, Proposition 4.2.3]. 

Here the Xi, Uj are formal variables, but we will assume that they are complex numbers 
in D when we apply formulas contained in this subsection. For each i — 1,2, ... ,n, we 
define qi{x) ^ qi{xi, . . . , x„) by 

= n (3-2) 

' l<k<n 

Using (3.1), we have 

qi{x)q2{x) ■ --qnix) = (-l)"("-i)/Met (-^ Y . (3.3) 

\ -^i-^i / i,j=i 

Recall the definition of Hafnians, which are sign-less analogs of Pfaffians. For a 2n x 2n 
symmetric matrix A — {aij)'^j^i, the Hafnian of A is defined by 

HiA^ ^2 ^r?{l)»?(2) «»7(3)J7(4) • • • ar?(2n-l)»7(2n) , (3.4) 
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summed over all permutations r) on {1,2, ...,2n} satisfying r]{2i — 1) < 'r]{2i) {i 
1, 2, . . . , n) and r]{l) < 77(8) < ■ ■ ■ < r]{2n - 1). 

A Pfaffian version of Cauchy's determinant identity is Schur's Pfaffian identity: 



2n 

I T.; T.- \ 

Pf 



n e 



» J / i,j=l l<i<j<2n * 

The following formula due to Ishikawa-Kawamuko-Okada [10] will be an important factor 
in our proofs of theorems. 

/ -r \ 2n / \ 2n 

lAj 7 lAj f -w- -w- I 1 . _ / <AJ •} iXj f 



l<i<j<2n f ^f ^i^i/i,j=l V(l ^i^j) J i,j=l 



3.2 Wick formula 

We recall the method for computations of expectations of polynomials in real Gaussian 
random variables. Let Yi, . . . , y„ be linear combinations of centered real Gaussian random 
variables. Then E[Yi • • • F„] = if n is odd, and 

E[Y,Y, •••!;] = m{E[Y,Y,m,^, (3.6) 

if n is even. For example, E[Y^Y2YsYi] ^ E[YiY2]E[Y-iYi]+E[YiYs\E[Y2Yi]+E[YiYi]E[Y2Y3]. 
See, e.g., survey [24] for details. 

4 Derivations of sign moments 

In this section, we provide a preliminary version of Pfaffian formulas for Theorem 5. 



4.1 Derivations for real Gaussian processes 

A derivation for sign moments of real Gaussian processes is given by a conditional expec- 
tation in the following way. 

Lemma 4.1. Let {X(t)}„i<j<i be a smooth real Gaussian process with covariance kernel 
K. Let {ti,t2, . . .,tn, and suppose that det K{i) — det(i^(ti, tj))'^^^^ 

does not vanish. Then, 

Qn 

E[sguX{ti) ■ • •sgnX(i„)sgnX(si) • • •sgnX(s^)] 



dtidt2---dt 

n/2 



[-\ E[X\t^) ■ --ritn) sgnX(si) • • -sgnXism) \ X{ti) = • • • = X{tn) = 0] 
X (detir(t))-^/2. 
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Proof. We will give a heuristic proof. The derivation of sgni is ^sgni = 25o(t), where 
5o(^) is Dirac's delta function at 0. Hence, if we abbreviate as Y{s) — sgnX(si) • • •sgnX(s^), 
then 

^^^^|l^E[sgnX(tO---sgnX(t„).y(^ 

=2"E[5o(X(ti))X'(ti) • ■ ■5o(X(t„))X'(y • Y{s)] 

=2"£;[X'(ti) • • • Y{s) I = • • • = = 0] •pt(0), 

where Pt(0) is the density of the Gaussian vector {X{ti), . . . , X{tn)) at (0, . . . , 0). Since 
pt(0) = (27r)-'^/2(^detir(t))-^/^ the claim follows. 

The above formal computation can be justified by using Watanabc's generalized Wiener 
functionals in the framework of Malliavin calculus over abstract Wiener spaces [22, 23]. □ 



4.2 Conditional expectations 

Recall the Gaussian power series / defined in (1.1). The process {/(t)}-i<t<i is centered 
real Gaussian with Cauchy covariance kernel 



/ ^ 1 



St 



The following identity in law is a crucial property which the Gaussian process with 
Cauchy covariance kernel enjoys. Set 

, , s — t 
^{s,t) 



1-st 

Lemma 4.2. For given {ti, t2, ■ ■ ■ , tn) £ X^, we have 

{f\f{ti)^---^f{tn)^0)^K;t)f (4.1) 
where /x(s,t) = IliLi /^l"^' ^0- Moreover, 

if, f{ti),i = 1, . . . , n\f{t,) = • • • = f{tn) = 0) ^ (//(•, t)f, qi{t)f{ti),i = 1, . . . , n), 
where qiit) = qi{ti, . . . ,tn) is defined in (3.2). 

Proof. If a Gaussian process X has a covariance kernel K[x, y), then that of the Gaussian 
process {X\X{t) = 0), i.e. X given X{t) = 0, is equal to K{x,y) - K{x,t)K{t,y) / K{t,t) 
whenever K{t, t) > 0. In the case of Cauchy kernel, we see that 

(^{x, y) - ^l^^l' = li{x, t)n{y, t)a{x, y). 
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This implies that = 0) = //(•, t)f as a process. Hence we obtain (4.1) by induction. 

As /' is a hnear functional of /, we also have the identity in law as a Gaussian system 

(/, /' I fih) = ■■■ = f{Q = 0) = t)f, /.'(•, t)f + /.(•, t)f). 

Since f^{ti,t) — and ii'{ti,t) — qi{t) for every i — l,2,...,n, we obtain the second 
equality in law. □ 

4.3 PfafRan expressions for derivations of signs 

The following lemma is a conclusion of Lemmas 4.1 and 4.2. 

Lemma 4.3. Let {t, s) = {ti, ^2, • • • , tn, Si, S2, ■ ■ ■ , Sm) ^ Xn+m- Then 

Qn 



dtidt2 ■■■dtn 



E[sgnf{ti) ■ ■ ■ sgn f{tn) sgn /(si) • • • sgn /(s^)] 



i=i j=i 

n/2 

(det • • ■ f{Q sgn/(si) ■ ■ ■ sgn /(s^)], 



X I - 



where ^{t) = {a{ti,tj))l^^,. 

Proof. From Lemma 4.2 and (3.3) we have 

E[f'{h) ■ ■ ■ f'itn) sgn /(si) ■ ■ ■ sgn /(s^) | f{h) = ■■■ = f{tn) = 0] 

m n 

j=i i=i 

n m 

= n n - • (-1)"^""'^/' det E(t) • £;[/(ii) • • • /(i„) sgn f{s,) • • • sgn /(s^)] . 

We have finished the proof by Lemma 4.1 with X{t) — f{t). □ 
Proposition 4.4. For t = (ti, . . . , t2n) G ^2n; we have 

Q2n 

-— -— ^[sgn/(ti)sgn/(t2) • • •sgn/(t2n)] 

atiat2 ■ ■ ■ Ot2n 
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Proof. Lemma 4.3 with m = and with the replacement n by 2n gives 

p)2n / 2 \ " 

^i?[sgn/(tO • • •sgn/(U] = (-1)'^ j {deimf"E[f{t,) ■ ■ ■ /(t^n)]- 
Here the Wick formula (3.6) gives 

Eim) ■ ■ ■ fit2n)] = mE[m)fitMu = Hf ((i - ut,)-%^, 

and Cauchy's determinant identity (3.1) gives 

(-l)'^(detE(t))V^= n 'Sn{t,-t^)-ll^==^- n 

l<i<i<2n i=l V '■i l<j<i<2n * ^ 



Hence we obtain 

g2n 



-E[sgn f{ti) sgn /(ti) • • • sgn f{t2n)] 



dtidt2 ■ ■ ■ dt2n 

^ ^ l<i<j<2n i=l V S l<i<j<2n * ^ * ^ «.:/=l 

The desired Pfaffian expression follows from the Pfafhan-Hafnian identity (3.5). □ 



5 Proof of Theorems 4 and 5 

5.1 Proof of Theorem 5 

Proposition 4.4 can be expressed as 

Q2n 



dtidt2 ■ ■ ■ dt2n 
2 



E[sgn f{ti) sgn /(ta) • • • sgn /(t 



2n) 



^ n sgn(t,-t,)-Pf(Kn(t„t,))2=, 

l<i<j<2n 

/ o \ " a2n 

(-) n «gnfe-*.).5j^^Pf(K.(,, . 



l<i<j<2n 



If we remove the differential symbol ^^^^ in the above equation, then we get the equality 
on Theorem 5. The goal of the present subsection is to prove that this observation is 
veritably true. 
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For each subset set I of {1,2,..., 2n}, we define the 2n x 2n skew symmetric matrix 
L''^ = lJ{t), the (i, j)-entry of which is 



Li 



^K22{ti,tj) iiielandjeP, 
K2i(ti, tj) = i-K22{ti, tj) ifieP and j e /, 

K22{U,tj) ifi,je/^. 

In particular, we put LI*^] = LMf 7 = {1, 2, . . . , A;} and = lJ>. 
Lemma 5.1. The following two claims hold true. 

1. For each A; = 0, 1, . . . , 2n - 1, Pf = Pf Ll'^+i] . 

2. Pf Ll'^l is skew symmetric in tfe+i, ■ ■ ■ , ^2n- 
Proof. Recall the definition of the Pfafhan 



(5.1) 



d 



d 



dtk+i 

For each i < j,we see that: 
if J = A; + 1, then 
if i = A; + 1, then 
if i,j 7^ A; + 1, then 
Hence we have 



PfLW = 5^(sgn77)-^J]L 



[k] 

■"7?(2i-l)r,(2i)- 



i=l 



d 



-L 



[k] 



d 



d 



L 



[k\ 



d 



_d 
W 



^22{tk+iitj) — 'K.i2{tk+i,tj) — 



fe+i 



-Lg = 0. 



d 



dt 



k+l 



Pf LW = ^(sgn,) niLKii,,,.., = Pf 



which is the first claim. 

Pfaffians are skew symmetric with respect to the change of the order of rows/columns, 
i.e., Pf (a^(j)^(j)) = (sgn?7)Pf A for any 2n x 2n skew symmetric matrix A = (aij) and a 
permutation rj on {1, 2, ... , 2n}. Hence the second claim follows from the definition of 
LW. □ 



Put 



■^2n 



{{ti, . . . , t2n) G ^2n | ^1 < • • • < t2n}- 
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Lemma 5.2. 



lim £;[sgn/(ii)---sgn/(i2n)] = 0. 
(ti,...,t2„)e^< 



Proof. If we put X{t) = VT^f{t), then 

^[sgn/(ti) • • ■sgn/(t2n)] = E[sgnX{ti) ■ • •sgnX(t2n)] 

and E[X{ti)X{tj)] = c(tj,tj) = ^ ' ■ Furthermore, since \imt^^^ic{ti,t2n) = (^i,2n, 

the random variable X{t2n) converges in distribution to a standard Gaussian variable 
independent of other Xi {i < 2n), which means that E[sgn X {ti) ■ ■ ■sgnX{t2n)] 0- □ 

Lemma 5.3. For each k — 0,1,2 . . . ,2n — 1, 

lim PfLl'^1^0. 

(tl, — ,t2n)G^2ri 

Proof. Taking the limit t2n ^ 1, each entry in the last row and column of L^^l converges 
to zero, and so does Pf Lt*^'. □ 

Lemma 5.4. Let {ti, . . . , t2n) G ^2n- ^^''^ ^^^^ k = 0,1, . . . , 2n, 

Qk 



E[sgn/(tO • ■ •sgn/(t2n)] = Pf LW- (5.2) 



dti--- dtk 

Proof. Consider the function Z^^^ on X2n defined by 

Qk 

Z^^\t,, . ..,t2n) ^^^-—^E[sgnf{t,) ■ ■■sgnf{t2n) 



-(l) n sgnfe-t,)-PfL['=]. 



l<i<j<2n 

Since Lt^"! = (Kii(ti, t^)) J^^, Proposition 4.4 implies that Z'^"! = on X2n- Let k G 
{0, 1, . . . , 2n — 1} and suppose that Zl'^+^l = on X2n- Our goal is to prove Z^''^ = on 

Prom the first statement of Lemma 5.1, -Q^^Z^'^^{ti, . . . , t2n) — Z^^^^\t-i, . .. , t2n), and 

hence our assumption implies Q^^Z^^\ti, . . . ,t2n) — 0. Therefore Z^^^ is independent of 

tk+i- Prom the second statement of Lemma 5.1, Z^^^ is symmetric in tk+i, ■ ■ ■ ,t2m and 
therefore Z^''^ is also independent of t2n- However, 

lim Zl^^l = 

t2n^l 
(tl,--->*2n)eX2„ 

by Lemmas 5.2 and 5.3. Hence, Z^''^ must be identically zero on X2„- □ 
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Proof of Theorem 5. Lemma 5.4 for A; = implies 

E[sgnf{h) ■ ■■sgnf{t2n)] = 0) Pf (K22(^i, 

for ti < • • • < t2n- Since Pf (1^22(^1, ^j))ij=i is skew symmetric in ti, . . . , t2n, 

n sgn(t,-i,)-Pf(K22(t„t,))g=i 

l<i<j<2n 

is symmetric and coincides with ^'[sgn f{ti) ■ ■ ■ sgn /(t2n)] on 3£2n- Thus we have obtained 
Theorem 5. □ 

The following corollary is a conclusion of Theorem 5. 

Corollary 7. For {ti, . . . ,t2n) G X2n o,nd a subset I = {ii < i2 < ■ ■ ■ < ik} in 
{l,2,...,2n}, 

£;[sgn f(h) • • • sgn f{t2n)] = f - I H " ' 



l<i<j<2n 



where L,^ is defined in (5.1). 

Proof Observe that g^-^^ Pf = Pf □ 

5.2 Proof of Theorem 4 
Lemma 5.5. For (ti, . . . , t„) G 

lim ■ • ■ lim 7- —E[sgnf{ti) ■ ■ •sgn/(i„) sgn/(si) • • •sgn/(s„)] 

si-s>ti+0 s„^tn+0 all ■ ■ ■ Otn 



Proof. Take — 1 < ti < si < t2 < S2 < ■ ■ ■ < < Sn < 1- Corollary 7 with / 
{1, 3, 5, . . . , 2n — 1} and with the replacement (ti, . . . , t2n) by (ti, Si, . . . , -Sn) gives 

-£^[sgn /(ii) • • • sgn f{tn) sgn /(si) • • • sgn 



dti--- dt 

TT J \K2i{Si,tj) K22{Si,Sj) 



Taking the limit Si — )■ ti,...,s„ — )■ t„, it converges to (|)" Pf (K(tj, tj))"^.^]^ for ti < 
■■■ < tn- From the symmetry for the achieved result holds true for every 

(ti, . . . , t„) e Xn- □ 



Correlation functions for zeros and PfafRans 



15 



Proof of Theorem 4- We use Lemma 4.3 with m — n. The identity in the lemma holds 
true iov —1 < ti < si < t2 < S2 < ■ ■ ■ < tn < Sn < i- Note that Y[i=i lYj=i ^S^l'^j ~ ^i) — 
(_l)n(n-i)/2_ Taking the limit ,si -^ti,...,Sn^ tn, 

Qri 

lim ••• lim — -— ^[sgn/(ti) ■ ■ ■sgn/(t„) sgn/(si) ■ ■ ■sgn/(s„)] 

si-s>ti+0 Sn^t„+0 Oil • • • <J'tn 

{<ieii:{t)f"E[\f{ti)---m\] 

for — 1 < ti < ■ ■ ■ < t„ < 1. From the symmetry for ti, . . . , t„, the above equation holds 
true for every (ti, . . . , t„) G X^- Combining this fact with Lemma 5.5, we obtain Theorem 
4. □ 



6 Proofs of Theorem 1, Corollary 2 and Corollary 3 
6.1 Proof of Theorem 1 

Hammersley's formula describes correlation functions of zeros of random polynomials, 
which was observed by Hammersley [8] and it is extended to Gaussian analytic functions 
as Corollary 3.4.2 in [9]. The following lemma is a real version of Hammersley's formula 
for correlation functions of Gaussian analytic functions. 

Lemma 6.1. Let X{t) be a random power series with independent real Gaussian coeffi- 
cients defined on an interval (—1, 1) with covariance kernel K . //det K{t) = det{K{ti, tj))^j^i 
does not vanish anywhere on X„, then the n-point correlation function for real zeros of f 
exists and is given by 

(, E[\X'{ti)---X'{tr:)\ I X(ti) = ■ ■ ■ = = 0] 

for t = {ti, . . . ,t„) e Xn- 

Proof. This can be proved in almost the same way as in the proof of (3.4.1) in Corollary 
3.4.2 in [9]. The only difference is that the exponent of \X'{ti) ■ ■ ■X'(t„)| is 1 in the case 
of real Gaussian coefficients instead of 2 in the complex case. This is due to the fact that 
the Jacobian determinant of F{t) = {X{ti), . . . , X{tn)) is equal to \X'{ti) ■ ■ when 
X is a real-valued differentiable function while \X'{ti) • • ■X'{tn)\'^ when X is complex- 
valued. □ 

Proof of Theorem 1. Prom Lemma 4.2 and (3.3) we have 

E[\f'{ti)---f'{t^)\ |/(tO = --- = /(tn)=0] 

= \qi{t)---q^{t)\E[\f{ti)---f{t^)\] 
= detE(t).£;[|/(tO---/(in)l], 
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and it follows from Lemma 6.1 that 

. . . = (27r)-"/2(det E(t))V2E[|/(ti) • • -/(^ |]. 
Hence Theorem 1 follows from Theorem 4. □ 

6.2 Proof of Corollaries 2 and 3 

Proof of Corollary 2. From (2.4) we observe that 

t) — \-\- |//(s, t)\c{s, t) arcsinc(s, t) — c{s, t^, 
and so R{s, s) — and R{s, ±1) = 1. A simple calculation yields 



9 . / ,N sgn(t - s 
— arcsmcls,!) = — 



and hence we obtain 



d s£fn( t — 9 ) 



where 

5f(s, t) := {c{s, tY arcsin c(s, t) — |/i.(s, arcsin c(s, t) + c(s, t)|/^(s, 
Since ±1) = and 



d . , — 4sgn(t — s) / M / NO • / N 
Q^dis^t) = ^— |/x(s,t)|c(s,i) arcsmc(s,t), 

we have g{s,t) > 0. This implies the claim. □ 



Proof of Corollary 3. The first equality immediately follows from ENr — f_^Pi{s)ds. 
Recall that 

/r r-r pr / f-r \ 2 

/ p2{s,t)dsdt + / pi{s)ds — I / pi{s)ds] . 
-r J —r J —T \J —r / 

We recall 2-correlation function 

P2(s, t) = 7r-'{Ki2(s, s)Ki2(t, t) - Ku{s, t)K22{s, t) + Ki2(s, t)K2i{s, t)} 
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as in (2.4). Taking the discontinuity of K22(s, t) sA, s — t into account and using integration 
by parts together with (2.3), we have 

/r 
Kll{s,t)K22{s,t)dt 
-r 

= l^ii. *)K22(s, + f)K22(s, *)l;..+o - /' Ki2(s, t)K,,(s, t)dt 

/r 
K,2is,t)K2,is,t)dt. 



It is easy to see that 



/r 
ds{Ki2{s, r)K22(s, r) - Ki2(s, -r)K22(s, -r)} 
-r 



i[K22(s,r)2 - K22(s, -rfY-, = K22{r,rf - K22(r, -r)' = 0(1). 



Hence, we obtain 



TT / Ki2{s,s)ds+ / Ki2(s,i)K2i(s,i)ciscin +0(1) 
= 27r- (^,rlogi±^-21ogl±^)+0(l). 

This imphes the assertion. □ 



7 Proof of Theorem 6 

7.1 Complex- valued Gaussian processes 

In this section, we assume that a process X = {X{X)}x^a is centered, i.e., E[X{X)] = 
for cadi A G A. Let X — {X{X)}x^a be a complex- valued Gaussian process in the sense 

that the real and imaginary parts form a real Gaussian process. We say that a complex- 
valued Gaussian process is a complex Gaussian process if the real and imaginary parts are 
mutually independent and have the same variance. 
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For a complex- valued Gaussian process X, we use three 2x2 matrices 



M{X,ii) - I 




M{X,n) = l 




M(A,/^) = I 





(E[ X{\)X{^,) ] E[X(A)X(/.)]\ 
\E[XiX)Xi^^)] E[XiX)Xi,^]J ' 

/E[X(A)X(/x)] E[ X(A)X(/x) ]\ 
U[^(A)X(/.)] £;[X(A)X(/.)]y' ' 

fE[^X{X)^X{ij,)] E[^X{X)^X{fj,)]\ 
\E[QX{X)?liX{n)] E[^X{X)^X{fj,)]) ■ 

For Ai, A2, . . . , A„ e A, the matrix (M(Ai, Xj))^j^^ is Hermitian, (M(Ai, Aj))"^-^;^ is complex 
symmetric, and (M(Ai, Xj))^^^^ is real symmetric. The real Gaussian vector 

(5RX(Ai), »X(Ai), . . . , 3fiX(A„), »^(A„)) 

has the covariance matrix (M(Aj, Aj))"^-^]^. It is easy to see that 

M(A,/.) = ^UM{\i^)U% U = (\-^ . (7.1) 



A (centered) complex-valued Gaussian process is uniquely determined by M or M. 
Lemma 7.1. For Ai, . . . , A„ G A, 

E[\x{x,)---x{xM-mn\,^jmj=i- (7.2) 

Proof. Let y„(A) = ?ftX(X) + a^X(X) for a e M. It follows from the Wick formula (3.6) 
that 

nuxM . . . UK)MK)] = Hf (^[;;{,';]^:{,';]j ^^^i^S^'^^^^ ■ (^-a) 

By analytic continuation, the formula (7.3) still holds for a, 6 e C. Therefore, by setting 
a — —b — we obtain the result. □ 

7.2 Conditional expectations for complex cases 

Throughout this section, we use the following notation. Put 

z-w 



1 — zw 

Let D+ = e D I > 0} and zi,...,Zn G D+. Moreover, we set as zj+n '■— 'z], 
j = 1, 2, . . . , n, and define qi{z) := qi{zx, Z2n) by (3.2) for i = 1, 2, . . . , 2n: 



■■^=Thf n ^ 



l<fe<2n ' 
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Let / be the Gaussian power series defined by (1.1). Then {/(-z)}zgD+ is a complex- 
valued Gaussian process with 



1 1 



\ l—zw 1—zw 

Lemma 7.2. For rj e D_|_, 

(/ I /(r/) = 0) = M-,^)M-,^)/- (7.4) 

Moreover, 

{f{z,),i = 1, 2, . . . , n I f{z,) = ■■■ = /(^J = 0) = iq^iz)f{z,),^ = 1, 2, . . . , n). (7.5) 

Proof. The real Gaussian vector {^f{z),Qf{w)) given ^/{rj) — Qf{rj) — has the co- 
variance matrix 

M/(^, w) - Mf{z, r))Mf{r), r))-^Mf{r), w) 
= ^U[Mf{z, w) - Mf{z, ri)Mf{ri, ri)-^Mf{ri, w)]U* 

by (7.1). A direct computation gives 

Mf{z,w) -Uf{z,rj)Mf{rj,rj)-^Mf{rj,w) = Uy^{z,w) 

with Yr^{z) = fj,{z,ri)fj,{z,fj)f{z), and we obtain (7.4). The remaining statement follows 
from (7.4) in a manner similar to the proof of Lemma 4.2. □ 

7.3 Correlation functions for complex zeros 

We finally compute the correlation function p^izi, . . . ,Zn) for complex zeros of /. Our 
starting point is the following Hammersley's formula (complex version), see [9] and com- 
pare with Lemma 6.1: 

. E[\f{z,)---nZr.)\'\f{z^) = ... = f{z^)=0] 
p^[Zi, ...,Zn)^ / ■ (7.Dj 

{27r)-^det{Mf{zi,Zj))l^^, 
Note that {27r)^"'[det(Mf{zi,Zj))]~^^^ is the density of the real Gaussian vector 

m(z,),^f(z,),...,m^n),^f(Zn)) 

at (0,0, ...,0,0). 
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Proposition 7.3. Let 

ThcTi 

(-1)" det M{z) ■ Hf M{z) 

P„{Zi, ...,Zn)- — ^====== . 

Proof. Let us compute the numerator on (7.6). Equation (7.5) gives 
E[\f{z,) ■ • •/'(^„)r I f{zi) = ■■■ = f{zn) = 0] = \q,{z) ■ ■■qr,{z)\'E[\f{z,) ■ --fizM- 



Here, since qj{z) — 5j+„(z) for j = 1, 2, . . . , n, it follows from (3.3) that 



2n 



\q,{z)---qn{z)\' = l[q^{z) = (-1)'^ det M(z). 



i=l 



Furthermore, from (7.2) we have 

E[\fiz,) . --fizM = m{Mf{z„Zj))l,^, = Hf M(z). 
On the other hand, the denominator on (7.6) is computed by using (7.1): 

det{M f{z,,Zj))l^^, = A-- det M{z). 
Consequently, we obtain the result from (7.6). □ 
Proof of Theorem 6. By the Cauchy determinant formula (3.1), 



detM(z)=n^ n (s) . 

i=l ^ l<i<j<2n ^ ^'^^ ^ 

In ^ 

1=1 V \Zi'\ l<i<i< 



Zi Zj 



l<i<j<2n 

By noting that Zi+n = (i = 1, 2, . . . , n), it is easy to see that 



1 Zj^Zj 



dot M{z) ^ -|^Nn(n-l)/2 TJ ^ "TT ^» ~ ^» ( TT ~ 



Since = ^f^-i for 'isz > 0, from Proposition 7.3 and Pfaffian-Hafnian identity (3.5) 
we see that 



(_l)"("-l)/2 " I / 2:- — z- \ 



2n 
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By changing rows and columns, we finally obtain 



1 



n 



1 



Pn{zi, ■ ■ ■ ,Zn) 



n 



n 



5j-Pf(K°(zi,z,.)) 




i=l 



□ 
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